
function phi = vercheil2(vec,l,dbete,t)
phi  = vec*cos(-vec*log(cos(deg2rad(dbete*t)))/(l*deg2rad(dbete)));
end